An Apparatus Having an Array of Piezoelectric Film Sensors 



for Measuring Parameters of a Process Flow Within a Pipe 



Cross-Reference To Related Patent Applications 
5 This application claims the benefit of U.S. Provisional Application No. 60/425,436 

(Cidra Docket No. CC-0538), filed November 12, 2002; and U.S. Provisional Application 
No. 60/426,724 (Cidra Docket No. CC-0554), filed November 15, 2002, all of which are 
incorporated herein by reference in their entirety. 

10 Technical Field 

This invention relates to an apparatus for measuring the parameters of a single phase 
and/or multiphase flow, and more particularly to an apparatus having an array of 
piezoelectric film sensors clamped or mounted onto a process flow pipe for measuring the 
speed of sound and/or vortical disturbances propagating in a single phase and/or multiphase 
15 flow to determine parameters, such as mixture quality, particle size, vapor/mass ratio, 

liquid/vapor ratio, mass flow rate, enthalpy and volumetric flow rate of the flow in the pipe, 
for example, by measuring acoustic and/or dynamic pressures. 

Background Art 

20 Numerous technologies have been implemented to measure volumetric and mass 

flow rates of fluids in industrial processes. Some of the more common approaches are 
based upon ultrasonic time of flight and/or Doppler effects, Coriolis effects, rotating 
wheels, electromagnetic induction, and pressure differentials. Each of these techniques has 
certain drawbacks. For example, invasive techniques that rely on insertion of a probe into 

25 the flow, or geometry changes in the pipe, may be disruptive to the process and prone to 
clogging. Other methods such as ultrasonics may be susceptible to air or stratified flow. 
Meters that use rotating wheels or moving parts are subject to reliability issues. Coriolis 
meters are limited when pipe diameters become large due to the increase in force required 
to vibrate the pipe. 

30 One such process fluid is a saturated vapor/liquid fluid mixture (e.g., steam). It 

would be advantageous to be able to measure the vapor quality of this fluid mixture. Vapor 
quality of a saturated vapor / liquid mixture is defined as ratio of the mass of the vapor 
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phase to the total mass of the mixture. Saturated mixtures exist at temperatures and 
pressures at which Hquid and vapor phases coexist. The temperatures and pressures at 
which the Hquid and vapor phases coexist lie under the "vapor bubble" on a phase diagram. 
The collection of points knovra as the saturated liquid line and the collections of points 
5 known as the saturated vapor line define the vapor bubble. These two lines connect at, what 
is termed, the critical point. Saturated mixtures exist only under the vapor bubble. For 
pressures and temperatures outside of the vapor bubble, the fluid exists as a single phase 
and the properties of that fluid, such as density, enthalpy, internal energy, etc., are uniquely 
defined by the pressure and temperature. For common fluids, such as water, these 

10 properties are tabulated as functions of pressure and temperatures and are available through 
a variety of references including a website hosted by NIST (ref: 
http://webbook.nist.gov/chemistrv/fluid/) . 

For fluids at pressures and temperatures that lie within the vapor bubble, the fluids 
represent mixtures of the liquid and vapor phase. Although the properties of both the vapor 

15 and liquid phases are well defined (and tabulated for known substances), the properties of 
the mixture are no longer uniquely defined as fimctions of pressure and temperature. In 
order to define the averaged properties of a saturated mixture, the ratio of the vapor and 
liquid components of the mixture must be defined. The quality of the mixture, in addition 
to the pressure and temperature, must be defined to uniquely determine the properties of the 

20 mixture. 

Measuring the average properties of a single or multi-phase process flow is 
important in many industrial application since it is the mass averaged properties of the 
working fluid that enter directly into monitoring the thermodynamic performance of many 
processes. For example, it is the difference in the flux of enthalpy of the steam mixture 

25 flowing into and exiting fi-om a turbine that determines the maximum mechanical work that 
can be extracted fi-om the working fluid, and thus is critical to determining component 
efficiency. However, if the steam entering or exiting the turbine were saturated, pressure 
and temperature measurement would not sufficient to determine the specific enthalpy, but 
rather, a measure of the quality of the steam would be required to uniquely define the 

30 thermodynamic properties of the saturated steam mixture. 
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Note that once the quality and pressure (or temperature) of a saturated mixture is 
defined, the thermodynamic properties of the mixture are defined through mixing laws 
provided the properties of the liquid and vapor sates are known. For example, measuring 
speed of sound enables one to determine quality, which in turn enables one to calculate 
5 enthalpy, density, and other properties of the mixture. In addition to measuring the specific 
enthalpy, a measurement of the total mass is also, in general, required to determine the flux 
of enthalpy. 

There are many other situations where knowing the quality of a saturated mixture is 
beneficial. For example, in a steam power plant, the quality of the steam within the steam 

10 turbine affects blade life. Generally it is desired to operate so the quality is as high as 
possible throughout the turbine to minimize liquid water drops that will erode the metal 
blades. Knowing the quality at the turbine inlet and exhaust (or at the exhaust only if the 
inlet is super-heated) provides a means to monitor the quality throughout the turbine. Also, 
to monitor plant performance so that it can be operated at optimum conditions and to 

15 identify degradation effects, the steam turbine thermal performance must be known. This 
requires the fluid enthalpy at the inlet and exhaust of each turbine to be known. If the fluid 
at either or both locations is saturated, pressure and temperature measurements alone will 
not be enough to determine the enthalpy. However if an additional measurement of quality 
is made the enthalpy is then defined. In addition, there may be other applications in 

20 refrigeration cycles. 

The ability to measure the flow rate and composition of the saturated vapor / liquid 
mixtures within the conduits is an important aspect of any system or strategy design to 
optimize the performance of a system based on saturated vapor / liquid mixtures. The 
industry recognizes this, and has been developing a wide variety of technologies to perform 

25 this measurement. These include probe based devices, sampling devices, Venturis and 
ultrasonic devices 

This invention provides an apparatus and method to measure homogeneous and/or 
non-homogeneous fluids used in industrial systems having various working fluids to 
determine various parameters of the process fluid, such as the volumetric flow of the fluid, 
30 the consistency or composition of the fluid, the density of the fluid, the Mach number of the 
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fluid, the size of particle flowing through the fluid, the air/mass ratio of the fluid and/or the 
percentage of entrained air/gas within a liquid or slurry. 

Here a novel approach to flow measurements is proposed which utilizes a non- 
intrusive, externally mounted sensing element that requires no moving parts and is highly 
5 reliable. This approach is based upon signal correlation and/or array processing techniques 
of unsteady pressure measurements induced in an array of externally mounted sensors. The 
piezo-film sensors clamped onto the outer surface of a pipe provides circumferential 
averaging of the unsteady pressures within the pipe and provide an inexpensive solution to 
accurately measuring the unsteady pressures. The piezo-film also have the advantage of 
10 being able to wrap around a substantial portion of the outer circumference of the pipe to 
provide circumferential averaging of the unsteady pressures with the pipe. 

Summary of the Invention 

Objects of the present invention include an apparatus having an array of 
15 piezoelectric film sensors mounted or clamped axially spaced to the outer surface of the 

pipe for measuring the imsteady pressures of a single and multi-phase process flows within 
a pipe to determine at least one parameter of the process flow. 

According to the present invention, an apparatus for measuring at least one 
parameter of a process flow flowing within a pipe. The apparatus includes at least two 
20 pressure sensors disposed on the outer surface of the pipe at different axial locations along 
the pipe. Each of the pressure sensors provides a respective pressure signal indicative of a 
pressure disturbance within the pipe at a corresponding axial position. Each of the pressure 
sensors includes a piezoelectric film sensor. A signal processor, responsive to said pressure 
signals, provides a signal indicative of at least one parameter of the process flow flowing 
25 within the pipe. 

The foregoing and other objects, features and advantages of the present invention 
will become more apparent in light of the following detailed description of exemplary 
embodiments thereof. 
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Brief Description of the Drawings 

Fig. la is a schematic illustration of an apparatus having an array of piezoelectric 
film sensors clamped onto the outer surface a pipe, in accordance with the present 
invention. 

5 Fig. lb is a schematic illustration of an apparatus having an array of piezoelectric 

film sensors mounted on the outer surface of a pipe, in accordance with the present 
invention. 

Fig. 2 is a cross-sectional view of a pipe and array of sensors showing the turbulent 
structures within the pipe, in accordance with the present invention. 
10 Fig. 3 is a cross-sectional view of a piezoelectric film sensor in accordance with the 

present invention. 

Fig. 4 is a top plan view of a piezoelectric film sensor in accordance with the present 
invention. 

Fig. 5 is a cross-sectional view of a portion of the piezoelectric film sensor and 
15 clamp, in accordance with the present invention. 

Fig, 6 is a cross-sectional view of a portion of the piezoelectric film sensor and 
clamp, in accordance with the present invention. 

Fig. 7 is a cross-sectional view of a portion of the piezoelectric film sensor and 
clamp, in accordance with the present invention. 
20 Fig. 8 is a side elevational view of a portion of the piezoelectric film sensor and 

clamp showing a step in the attachment of the clamp to the pipe, in accordance with the 
present invention. 

Fig. 9 is a side elevational view of a portion of the piezoelectric film sensor and 
clamp, in accordance with the present invention. 
25 Fig. 10 is a perspective view of a plurality of piezoelectric film sensors clamped to a 

pipe having covers disposed thereover, in accordance with the present invention. 

Fig. 1 1 is a cross sectional end view of a piezoelectric film sensor clamped to a pipe, 
in accordance with the present invention. 

Fig. 12 is a block diagram of a probe for measuring the speed of sound propagating 
30 through a process flow flowing within a pipe, in accordance with the present invention. 
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Fig. 1 3 is a plot showing the standard deviation of sound speed versus frequency for 
various arrays of process flow parameter measurement system, in accordance with the 
present invention. 

Fig. 14 is a kco plot of data processed from an array of pressure sensors use to 
measure the speed of sound propagating through a saturated vapor/liquid mixture flowing in 
a pipe, in accordance with the present invention. 

Fig. 15 is a block diagram of an apparatus for measuring the vortical field of a 
process flow within a pipe, in accordance with the present invention. 

Fig. 16 is a kcj plot of data processed from an apparatus embodying the present 
invention that illustrates slope of the convective ridge, and a plot of the optimization 
function of the convective ridge, in accordance with the present invention. 

Fig. 17 is a functional flow diagram of an apparatus embodymg the present 
invention. 

Best Mode for Cartying Out the Invention 

Referring to Fig. la, an apparatus, generally shown as 10, is provided to sense and 
determine specific characteristics or parameters of a single phase fluid 12 (e.g., gas and 
liquid) and/or a multi-phase mixture 12 (e.g., process flow) flowing through a pipe. The 
multi-phase mixture may be a two-phase liquid/vapor mixture, a solid/vapor mixture or a 
solid/liquid mixture, gas entrained liquid or even a three-phase mixture. As will be 
described in greater detail, the apparatus measures the speed of sound propagating through 
the fluid or multiphase mixture flow to determine any one of a plurality of parameters of the 
flow, such as the steam quality or * Vetness", vapor/mass ratio, liquid/solid ratio, the 
volumetric flow rate, the mass flow rate, the size of the suspended particles, density, gas 
volume fraction, and the enthalpy of the flow. Additionally, the apparatus 10 is capable of 
measuring the unsteady pressure disturbances (e.g., vortical effects, density changes) of the 
flow passing through the pipe to determine the velocity of the flow, and hence the 
volumetric flow rate of the flow. 

Fig. la illustrates a schematic drawing of the apparatus 10 that includes a sensing 
device 16 comprising an array of pressure sensors (or transducers) 18-21 spaced axially 
along the outer surface 22 of a pipe 14, having a process flow propagating therein. The 
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pressure sensors measure the unsteady pressures produced by acoustical and/or vortical 
disturbances within the pipe, which are indicative of a parameter of the single phase fluid or 
multiphase mixture 12. The output signals (Pi - P4) of the pressure sensors 18-21 are 
provided to a processing imit 24, which processes the pressure measurement data and 
determines at least one parameter of the flow. Such as, the characteristics and parameters 
determined may include the volumetric flow of the flow, the consistency or composition of 
the flow, the density of the mixture, the Mach number of the flow, the size of particle 
flowing through the mixture, the air/mass ratio of the mixture and/or the percentage of 
entrained air or gas within the mixture. 

In an embodiment of the present invention shown in Fig. la, the apparatus 10 has 
four pressure sensors 18-21 disposed axially along the pipe 14 for measuring the unsteady 
pressure P1-P4 of the fluid or mixture 12 flowing therethrough. The apparatus 10 has the 
ability to measure the volumetric flow rate and other flow parameters using one or both of 
the following techniques described herein below: 

1) Determining the speed of sound of acoustical disturbances or sound waves 
propagating through the flow 12 using the array of pressure sensors 18-21, 
and/or 

2) Determining the velocity of vortical disturbances or "eddies" propagating 
through the flow 12 using the array of pressure sensors 18-21. 

Generally, the first technique measures imsteady pressures created by acoustical 
disturbances propagating through the flow 12 to determine the speed of sound (SOS) 
propagating through the flow. Knowing the pressure and/or temperature of the flow and the 
speed of soimd of the acoustical disturbances, the processing unit 24 can determine the mass 
flow rate, the consistency of the mixture (i.e., the mass/air ratio, the mass/liquid ratio, the 
liquid/air ratio), the volumetric flow rate, the density of the mixture, the enthalpy of the 
mixture, the Mach number of the mixture, the size of the particles within a mixture, and 
other parameters, which will be described in greater detail hereinafter. 

The apparatus in Fig. la also contemplates providing one or more acoustic sources 
27 to enable the measurement of the speed of soimd propagating through the flow for 
instances of acoustically quiet flow. The acoustic sources may be disposed at the input end 
of output end of the probe, or at both ends as shown. One should appreciate that in most 
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instances the acoustics sources are not necessary and the apparatus passively detects the 
acoustic ridge provided in the flow 12. 

The second technique measures the velocities associated with unsteady flow fields 
and/or pressure disturbances created by vortical disturbances or "eddies" 1 1 8 to determine 
5 the velocity of the flow 12. The pressure sensors 18-21 measure the unsteady pressures Pi- 
P4 created by the vortical disturbances as these disturbances convect within the flow 12 
through the pipe 14 in a known manner, as shown in Fig. 2. Therefore, the velocity of these 
vortical disturbances is related to the velocity of the mixture and hence the volumetric flow 
rate may be determined, as will be described in greater detail hereinafter. 
10 In one embodiment of the present invention as shown in Figs, la, lb and 2, each of 

the pressure sensors 18-21 may include a piezoelectric film sensor 30 as shown in Figs. 3 
and 4 to measure the unsteady pressures of the mixture 12 using either technique described 
hereinbefore. 

As best shown in Figs. 3 and 4, the piezoelectric film sensors 30 include a 

15 piezoelectric material or film 32 to generate an electrical signal proportional to the degree 
that the material is mechanically deformed or stressed. The piezoelectric sensing element is 
typically conformed to allow complete or nearly complete circumferential measurement of 
induced strain to provide a circumferential-averaged pressure signal. The sensors can be 
formed fi-om PVDF films, co-polymer films, or flexible PZT sensors, similar to that 

20 described in "Piezo Film Sensors Technical Manual" provided by Measurement Specialties, 
Inc., which is incorporated herein by reference. A piezoelectric film sensor that may be 
used for the present invention is part number 1-1002405-0, LDT4-028K, manufactured by 
Measurement Specialties, Inc. 

Piezoelectric film ("piezofilm"), like piezoelectric material, is a dynamic material 

25 that develops an electrical charge proportional to a change in mechanical stress. 

Consequently, the piezoelectric material measures the strain induced within the pipe 14 due 
to unsteady pressure variations (e.g., vortical and/or acoustical) within the process mixture 
12. Strain within the pipe is transduced to an output voltage or current by the attached 
piezoelectric sensor. The piezoelectrical material or film may be formed of a polymer, such 

30 as polarized fluoropolymer, polyvinylidene fluoride (PVDF). 
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Figs. 3 and 4 illustrate a piezoelectric film sensor (similar to the sensor 18 of Fig. 
la), wherein the piezoelectric film 32 is disposed between and pair of conductive coatings 
34,35, such as silver ink. The piezoelectric film 32 and conductive coatings 34,35 are 
coated onto a protective sheet 36 (e.g., mylar) with a protective coating 38 disposed on the 
opposing side of the upper conductive coating. A pair of conductors 40,42 is attached to a 
respective conductive coating 34,35. 

The thickness of the piezoelectric film 32 may be in the range of 8 um to 
approximately 110 um. The thickness is dependent on the degree of sensitivity desired or 
needed to measure the unsteady pressures within the pipe 14. The sensitivity of the sensor 
30 increases as the thickness of the piezoelectric film increases. 

The advantages of this technique of clamping the PVDF sensor 30 onto the outer 
surface of the pipe 14 are the following: 

1. Non-intrusive flow rate measurements 

2. Low cost 

3. Measurement technique requires no excitation source. Ambient flow noise is 
used as a source. 

4. Flexible piezoelectric sensors can be mounted in a variety of configurations 
to enhance signal detection schemes. These configurations include a) co- 
located sensors, b) segmented sensors with opposing polarity configurations, 
c) wide sensors to enhance acoustic signal detection and minimize vortical 
noise detection, d) tailored sensor geometries to minimize sensitivity to tube 
modes, e) differencing of sensors to eliminate acoustic noise from vortical 
signals. 

5 . Higher Operating Temperatures ( 1 25C) (co-polymers) 

Referring to Fig. lb, the piezoelectric film sensors 30 may be mounted directly onto 
the outer diameter of the pipe 14 by epoxy, glue or other adhesive. 

Alternatively, as shown in Figs, 5-9, the piezoelectric film sensor 30 is be adhered 
or attached to a strap 72 which is then clamped (or strapped) onto the outer surface of the 
pipe 14 at each respective axial location, similar to that described in U.S. Provisional 
Application No. 60/425,436 (Cidra Docket No. CC-0538), filed November 12, 2002; and 
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U.S. Provisional Application No. 60/426,724 (Cidra Docket No. CC-0554), which are 
incorporated herein by reference. 

As shown in Fig 5, the piezoelectric film sensor 30 is attached to the outer surface 
73 of the strap in relation to the pipe 14. The conductive insulator 36 is attached to the 
outer surface of the strap by double side tape or any other appropriate adhesive. The 
adhesive is preferably flexible or compliant but minimizes creep between the strap and 
piezoelectric film sensor during the operation of the sensor 30. The length of the strap is 
substantially the same as the circumference of the pipe 14. The piezoelectric film sensor 
may extend over the substantial length of the strap or some portion less than the strap. In 
the embodiment shown in Fig.5, the piezoelectric film sensor 30 extends substantially the 
length of the strap 72 to provide a circumferentially averaged pressure signal to the 
processing unit 24. 

Referring to Figs. 6 and 7, an attachment assembly 75 comprising a first attachment 
block 76, a second attachment block 77 and a spacer 78 disposed therebetween, which are 
welded together to provide slots 79 between each of the attachment blocks and the spacer. 
The slots receive respective ends of the strap 72 to secure the ends of the strap together. One 
end of the strap 72 and the pair of conductors 40,42 are threaded through the slot disposed 
between the first attachment block and the spacer. The strap 72 and conductors 40,42 are 
secured to the attachment assembly by a pair of fasteners 80. The other end of the strap is 
threaded through the slot disposed between the spacer and the second attachment block. 
Referring to Fig. 8, the other end of the strap is pull tightly between the spacer and the 
second attachment to draw up and take-up the tension and securely clamp the strap to the 
pipe 14. As shown in Fig. 9, a set screw 81 within the second attachment block is tighten, 
which then pierces the other end of the strap to secure it to the attachment assembly. The 
excess portion of the other end of the strap is then cut off. The piezoelectric film sensor 
may then be covered with a copper sheet to provide a grounding shield for EMI or other 
electrical noise. 

While the piezoelectric film sensor 30 was mounted to the outer surface of the straps 
72, the present invention contemplates the piezoelectric film sensor may be mounted to the 
inner surface of the strap, thereby resulting in the piezoelectric sensor being disposed 
between the strap and the outer surface of the pipe 14. 
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Fig. 10 illustrates a protective cover 82, having two halves, clamped onto the pipe 
and secured together over each of the piezoelectric film sensors 30 and straps 72. The 
protective cover is formed from aluminum with thermal fins 83 molded therein to assist 
with dissipating heat away from the sensors. The cover fiirther includes an insulative 
portion 84 disposed between the pipe and the aluminum portion of the cover. The 
conductors and wiring thereto pass through a conduit 84 that extends between each of the 
covers. 

While the present invention illustrates separate covers for each sensor, the present 
invention contemplates a single cover that covers all the sensors. 

Referring to Figs. 13 and 14, an apparatus 1 10, similar to apparatus 10 of Fig. la, 
embodying the present invention is provided that measures at least one 
parameter/characteristic of a single and/or multiphase flow 12 flowing within a pipe 14. 
The apparatus may be configured and programmed to measure the speed of sound 
propagating through the flow 12 or measure the vortical disturbances propagating through 
the flow 12. In some instances, the apparatus 10 may be configured to measure both the 
speed of sound and the vortical disturbances. Depending on the configuration or 
embodiment, the apparatus can measure at least one of the following parameters of the flow 
12: the wetness or steam quality (volumetric phase fi:action), the volumetric flow rate, the 
size of the liquid particles, the mass flow, the enthalpy and the velocity of the mixture. To 
determine any one of these parameters, the apparatus 110 measures the unsteady pressures 
created by the speed of sound (SOS) and/or the vortical disturbances propagating through 
the single phase or multiphase flow 12 flowing in the pipe 14, which will be described in 
greater detail hereinafter. 

The type of unsteady pressure measurement being made determines the spacing of 
the sensors. Measurement of unsteady vortical pressures require sensors spacing less than 
the coherence length of the vortical disturbances which is typically on the order of a pipe 
diameter. Correlation or array processing of the unsteady vortical pressure measurements 
between sensors is used to determine the bulk flow rate of the process mixture, which will 
be described in greater detail hereinafter. 

Mass flow rates and other parameters are determined by measuring the speed of 
sound propagating within the process mixture 12. These parameters are determined by 
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correlating or array processing unsteady pressure variations created by acoustic disturbances 
within the process mixture. In this case, the wavelength of the measured acoustic signal 
determines the sensor spacing. The desired wavelength of the measured acoustic signal is 
dependent upon the dispersion of particles in the mixture flow, which is dependent on the 
5 particle size. The larger the particle size is the longer the sensing device of the aperture. 

As described hereinbefore, the apparatus 1 10 of the present invention may be 
configured and programmed to measure and process the detected unsteady pressures Pi(t) - 
PN(t) created by acoustic waves and/or vortical disturbances, respectively, propagating 
through the mixture to determine parameters of the mixture flow 12. One such apparatus 10 

10 is shown in Fig. 12 that measures the speed of sound (SOS) of one-dimensional sound 
waves propagating through the vapor/liquid mixture to determine the composition the 
mixture. The apparatus 1 10 is also capable of determining the average size of the droplets, 
velocity of the mixture, enthalpy, mass flow, steam quality or wetness, density, and the 
volumetric flow rate of the single or multi-phase flow 12. It is known that sound propagates 

15 through various mediums at various speeds in such fields as SONAR and RADAR fields. 
The speed of sound propagating through the flow 12 within the pipe 14 may be determined 
using a number of known techniques, such as those set forth in U.S. Patent Application 
Serial No. 09/344,094, entitled "Fluid Parameter Measurement in Pipes Using Acoustic 
Pressures^ filed June 25, 1999, now US 6,354,147; U.S. Patent Application Serial No. 

20 09/729,994, filed December 4, 2002, now US 6,609,069; and U.S. Patent Application Serial 
No. 10/007,749, entitled "Fluid Parameter Measurement in Pipes Using Acoustic 
Pressures", filed November 7, 2001, each of which are incorporated herein by reference. 

In accordance with the present invention, the speed of sound propagating through 
the process flow 12 is measured by passively listening to the flow with an array of imsteady 

25 pressure sensors to determine the speed at which one-dimensional compression waves 
propagate through the flow 12 contained within the pipe 14. 

As shown in Fig. 13, the apparatus 110 has an array of at least three acoustic 
pressure sensors 1 15,1 16,1 17, located at three locations xi,X2,X3 axially along the pipe 14. 
One will appreciate that the sensor array may include more than three pressure sensors as 

30 depicted by pressure sensor 1 18 at location xn- The pressure generated by the acoustic 

waves may be measured through pressure sensors 115-118. The pressure sensors 15-18 
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provide pressure time-varying signals Pi(t),P2(t),P3(t),PN(t) on lines 120,121,122,123 to a 
signal processing unit 130 to known Fast Fourier Transform (FFT) logics 126,127,128,129, 
respectively. The FFT logics 126 - 129 calculate the Fourier transform of the time-based 
input signals Pi(t) - PwCt) and provide complex frequency domain (or frequency based) 
5 signals Pi(co),P2(co),P3((o),Pn((o) on lines 132,133,134,135 indicative of the frequency 
content of the input signals. Instead of FFT's, any other technique for obtaining the 
frequency domain characteristics of the signals Pi(t) - PN(t), may be used. For example, the 
cross-spectral density and the power spectral density may be used to form a frequency 
domain transfer fiinctions (or frequency response or ratios) discussed hereinafter. 

10 The frequency signals Pi(a)) - Pn(co) are fed to amix-Mx Calculation Logic 138 which 

provides a signal to line 40 indicative of the speed of sound of the multiphase mixture amix 
(discussed more hereinafter). The amix signal is provided to map (or equation) logic 142, 
which converts amix to a percent composition of a mixture and provides a %Comp signal to 
line 44 indicative thereof (as discussed hereinafter). Also, if the Mach number Mx is not 

IS negligible and is desired, the calculation logic 138 may also provide a signal Mx to line 46 
indicative of the Mach number Mx, 

More specifically, for planar one-dimensional acoustic waves in a homogenous 
mixture, it is known that the acoustic pressure field P(x,t) at a location x along the pipe 14, 
where the wavelength X of the acoustic waves to be measured is long compared to the 

20 diameter d of the pipe 14 (i.e., 7Jd »1), may be expressed as a superposition of a right 
traveling wave and a left traveling wave, as follows: 



where A,B are the frequency-based complex amplitudes of the right and left traveling 
waves, respectively, x is the pressure measurement location along a pipe 14, co is frequency 
25 (in rad/sec, where co=27if), and kr,ki are wave numbers for the right and left traveling waves, 
respectively, which are defined as: 




Eq. 1 




Eq.2 



CiDRA Corporation Proprietary 
CC-0676 



- 13- 



EV 286928210 US 



where amix is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), and 
Mx is the axial Mach number of the flow of the mixture within the pipe, where: 

M^^^ Eq.3 

where Vmix is the axial velocity of the mixture. For non-homogenous mixtures, the axial 
5 Mach number represents the average velocity of the mixture and the low frequency acoustic 
field description remains substantially unaltered. 

The data from the array of sensors 115-118 may be processed in any domain, 
including the frequency/spatial domain, the temporal/spatial domain, the temporal/wave- 
number domain or the wave-number/frequency (k-co) domain. As such, any known array 
10 processing technique in any of these or other related domains may be used if desired, 
similar to the techniques used in the fields of SONAR and RADAR. 

Also, some or all of the functions within the signal processing unit 130 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 
implemented using analog and/or digital hardware, having sufficient memory, interfaces, 
15 and capacity to perform the fiuictions described herein. 

Acoustic pressure sensors 1 15 - 118 sense acoustic pressure signals that, as 
measured, are lower frequency (and longer wavelength) signals than those used for 
ultrasonic probes of the prior art, and thus the current invention is more tolerant to 
inhomogeneities in the flow, such as time and space domain inhomogeneities within the 
20 flow. 

It is within the scope of the present invention that the pressure sensor spacing may 
be known or arbitrary and that as few as two sensors are required if certain information is 
known about the acoustic properties of the process flow 12. The pressure sensors are 
spaced sufficiently such that the entire length of the array (aperture) is at least a significant 
25 fraction of the measured wavelength of the acoustic waves being measured. The acoustic 
wavelength to be measured in a mixture is a fimction of at least the size and mass of the 
droplets/particles, and the viscosity of the vapor. The greater the size and mass of the 
droplets and/or the less viscous the vapor, the greater the spacing of the sensors is needed. 
Conversely, the smaller the size and mass of the droplets/particles and/or the more viscous 
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the vapor, the shorter the spacing of the sensors is needed. For single phase flow, the 
acoustic wavelength is a function of the type or characteristics of flow 12. 

Assuming that the droplets/particles of the mixture are small enough and the 
acoustic frequencies and the frequencies of perturbations associated with the acoustics are 
5 low enough for the droplets/particles of liquid to exhibit negligible slip (both steady and 
unsteady), the sound speed can be assumed to be substantially non-dispersive (that is 
constant with frequency) and the volumetric phase fraction of the mixture could be 
determined through the Wood equation: 

N 

10 Pn,ix = Z ^iPi 

_L_ = |;_^ 

Pmix^lix f=l 
N 

1=1 

15 

For one-dimensional waves propagating, the compliance introduced by the pipe ( in 
this case a circular tube of modulus E, radius R and wall thickness t) reduces the measured 
sound speed from the infinite dimensional soimd speed. The effect of the conduit is given 
by the following relationship: 

20 

= — +0" where o*s 

Pmix^measund Pmix^mix 

Utilizing the relations above for a vapor/liquid mixture, the speed at which sound 
travels within the representative vapor/liquid mixture is a function of vapor/liquid mass 
25 ratio. The effect of increasing liquid fraction, i.e. decreasing vapor/liquid ratio, is to 

decrease the sound speed. Physically, adding liquid droplets effectively mass loads the 
mixture, while not appreciably changing the compressibility of the air. Over the parameter 
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range of interest, the relation between mixture sound speed and vapor/liquid ratio is well 
behaved and monatomic. 

While the calibration curves based on predictions from first principles are 
encouraging, using empirical data mapping from sound speed to vapor/liquid ratio may 
5 result in improved accuracy of the present invention to measure the vapor/liquid fractions of 
the mixture. 

The sound speed increases with increasing frequency and asymptotes toward a 
constant value. The sound speed asymptote at higher frequency is essentially the sound 
speed of air only with no influence of the suspended liquid droplets. Also, it is apparent 

10 that the sound speed of the vapor/liquid mixture has not reached the quasi-steady limit at 
the lowest frequency for which sound speed was measured. The sound speed is continuing 
to decrease at the lower frequency limit. An important discovery of the present invention is 
that the speed at which sound propagates through droplets suspended in a continuous vapor 
is said to be dispersive. As defined herein, the speed at which acoustic waves propagate 

15 through dispersive mixtures varies with frequency. 

Measuring the sound speed of a saturated vapor/liquid mixture 12 at progressively 
lower and lower frequencies becomes inherently less accurate as the total length of the array 
of pressure sensors 1 15 - 118 (Axapemire), which define the aperture of the array, becomes 
small compared to the wavelength of the acoustics. In general, the aperture should be at 

20 least a significant fraction of a wavelength of the sound speed of interest. Consequently, 

longer arrays are used to resolve sound speeds at lower frequencies, which will be described 
m greater detail hereinafter. As shown in Fig. 14, the standard deviation associated with 
determining the speed of sound in air is shown as a fiinction of frequency for three arrays of 
varying aperture, namely 1.5 ft, 3 ft and 10ft. 

25 For accurately measuring sound speeds at ultra-low frequencies, the data suggests that 

utilizing a quasi-steady model to interpret the relationship between sound speed, measured 
at frequencies above those at which the quasi-steady model is applicable, and the liquid-to- 
vapor ratio would be problematic, and may, in fact, be impractical. Thus, the key to 
understanding and interpreting the composition of vapor/liquid mixtures through soimd 

30 speed measurements lies in the dispersive characteristics of the vapor/liquid mixture, which 
is described in greater detail in U.S. Patent Application Serial No. 10/412,839, filed April 
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10, 2003; U.S. Patent Application Serial No. 10/349,716, filed January 3, 2003; and U.S. 
Patent Application Serial No. 10/376,427, filed February 26, 2003, which are all 
incorporated herein by reference. 

In accordance with the present invention the dispersive nature of the system utilizes 
a first principles model of the interaction between the vapor and liquid droplets. This model 
is viewed as being representative of a class of models that seek to accoimt for dispersive 
effects. Other models could be used to account for dispersive effects without altering the 
intent of this disclosure (for example, see the paper titled "Viscous Attenuation of Acoustic 
Waves in Suspensions" by R.L. Gibson, Jr. and M.N. Toksoz), which is incorporated herein 
by reference. The model allows for slip between the local velocity of the continuous vapor 
phase and that of the droplets. The drag force on the droplets by the continuous vapor is 
modeled by a force proportional to the difference between the local vapor velocity and that 
of the liquid droplets and is balanced by inertial force: 



where K = proportionality constant, Uf = fluid velocity. Up = liquid droplet velocity, pp = 
liquid droplet density and Vp = particle volume. 

The effect of the force on the continuous vapor phase by the liquid droplets is modeled as a 
force term in the axial momentum equation. The axial momentum equation for a control 
volume of area A and length Ax is given by: 



where P = pressure at locations x and Ax, <t>p = volimie fraction of the liquid droplets, pf = 
vapor density. 

The droplet drag force is given by: 




dt 
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where Ca = drag coefficient, Ap = frontal area of liquid droplet and pf= vapor density. 

5 Using Stokes law for drag on a sphere at low Reynold's number gives the drag coefficient 
as: 



C = — = ^"^^ 
' Re p,[U,-U^)D^ 



10 where Dp = droplet diameter and \i = vapor viscosity. 



Solving for K in this model yields: 



15 



Using the above relations and 1 -dimensional acoustic modeling techniques, the following 
relation can be derived for the dispersive behavior of an idealized vapor/liquid mixture. 



20 



I + - 



Pf 



l + o 



25 



In the above relation, the fluid SOS, density (p) and viscosity (0) are those of the pure phase 
fluid, Vp is the volume of individual droplets and <t>p is the volumetric phase fraction of the 
droplets in the mixture. These relationships to determine droplet size and liquid to vapor 
mass ratio are described in U.S. Patent Application Serial No. 10/412,839, filed April 10, 
2003; U.S. Patent Application Serial No. 10/349,716, filed January 3, 2003; and U.S. Patent 
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Application Serial No. 10/376,427, filed February 26, 2003, which are all incorporated 
herein by reference. 

The apparatus 10 further includes the ability to measure of volumetric flow rate of 
the mixture by comparing the difference of the speed of one dimensional sound waves 
propagating with and against the mean flow. 

This method of determining the volumetric flow rate of the flow 12 relies on the 
interaction of the mean flow with the acoustic pressure field. The interaction results in 
sound waves propagating with the mean flow traveling at the speed of sound (if the 
vapor/liquid mixture were not flowing) plus the convection velocity and, conversely, sound 
waves traveling against the mean flow propagating at the speed of sound minus the 
convection velocity. That is, 

ClR = amix + u 

where aR = velocity of a right traveling acoustic wave relative to a stationary observer (i.e. 
the tube 14), at = velocity of a left traveling acoustic wave apparent to a stationary 
observer, amix == speed of sound traveling through the mixture (if the mixture was not 
flowing) and u = the mean flow velocity (assumed to be flowing fi-om left: to right in this 
instance). Combining these two equations yields an equation for the mean velocity. 




Therefore, by measuring the propagation velocity of acoustic waves in both directions 
relative to the pipe 14 as described hereinbefore, the mean flow velocity can be calculated 
by multiplying the mean flow velocity by the cross-sectional area of the pipe 14. 

The practicality of using this method to determine the mean flow is predicated on 
the ability to resolve the sound speed in both directions with sufflcient accuracy to 
determine the volumetric flow. 
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For the sound speed measurement, the apparatus 10 utilizes similar processing 
algorithms as those employed herein before, and described in greater detail hereinafter. The 
temporal and spatial frequency content of soimd propagating within the pipe 14 is related 
through a dispersion relationship. 

5 o)- 

The wave number is k, which is defined as k=27c/>-, co is the temporal frequency in rad/sec, 
and amix is the speed at which sound propagates within the process piping. For this cases 
where soimd propagates in both directions, the acoustic power is located along two acoustic 
ridges, one for the sound traveling with the flow at a speed of amix + Vmix and one for the 

10 sound traveling against the flow at a speed of amix - Vmix- 

Further, Fig. 14 illustrates the ability of the present invention to determine the 
velocity of a fluid moving in a pipe. The color contours represent the relative signal power 
at all combinations of frequency and wavenumber. The highest power "ridges" represent 
the acoustic wave with slope of the ridges equal to the propagation speed. Note that the 

1 5 acoustic ridges "wrap" to the opposite side of the plot at the spatial Nyquist wavenumber 
equal to ±3.14 in this case (i.e. the acoustic ridge that slopes up and to the right starting at 
the bottom of the plot, the right-side ridge, wraps to the left side of the plot at approximately 
550 Hz and continues sloping up and to the right). The dashed lines show the best-fit two- 
variable maximization of the power with the two variables being sound speed and flow 

20 velocity. The right-side ridge represents the acoustic wave traveling in the same direction 
as the bulk flow and therefore its slope is steeper than the left-side ridge that represents the 
acoustic wave traveling in the opposite direction of the flow. This indicates that the 
acoustic wave traveling in the same direction of the flow is traveling faster than the acoustic 
wave traveling in the opposite direction of the flow relative to the stationary sensors located 

25 on the probe. 

Referring to Fig. la, an apparatus 10 embodying the present invention includes the 
ability to measure volumetric flow rate of the mixture by measuring the unsteady pressures 
generated by vortical disturbance propagating in the mixture. The apparatus 10 uses one or 
both of the following techniques to determine the convection velocity of the vortical 
30 disturbances within the process flow 12 by: 
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1) Cross-correlating xinsteady pressure variations using an array of unsteady pressure 
sensors. 

2) Characterizing the convective ridge of the vortical disturbances using an array of 
unsteady pressure sensors. 

The overwhelming majority of industrial process flows involve turbulent flow. 
Turbulent fluctuations within the process flow govern many of the flow properties of 
practical interest including the pressure drop, heat transfer, and mixing. For engineering 
applications, considering only the time-averaged properties of turbulent flows is often 
sufficient for design purposes. For sonar flow metering technology, understanding the time- 
averaged velocity profile in turbulent flow provides a means to interpret the relationship 
between speed at which coherent structures convect and the volumetrically averaged flow 
rate. 

From the saturated vapor / liquid mixture mechanics perspective, this method relies 
on the ability of the apparatus 10 to isolate the convective pressure field (which convects at 
or near the mean velocity of the saturated vapor / liquid mixture) fi-om the acoustic pressure 
field (which propagates at the at the speed of sound). In this sense, the velocity 
measurement is independent of the sound speed measurement. 

For turbulent flows 12, the time-averaged axial velocity varies with radial position 
fi-om zero at the wall to a maximum at the centerline of the pipe 14. The flow near the wall 
is characterized by steep velocity gradients and transitions to relatively imiform core flow 
near the center of the pipe 14. Fig. 2 shows a representative schematic of a velocity profile 
and coherent vortical flow structures 188 present in fiilly developed turbulent flow 12. The 
vortical structures 188 are superimposed over time averaged velocity profile within the pipe 
14 and contain temporally and spatially random fluctuations with magnitudes typically less 
than 10% percent of the mean flow velocity. 

From a volumetric flow measurement perspective, the volumetrically averaged flow 
velocity is of interest. The volumetrically averaged flow velocity, defined as V=Q/A, is a 
usefiil, but arbitrarily defined property of the flow. Here, A is the cross sectional area of the 



CiDRA Corporation Proprietary 
CC-0676 



-21 - 



EV 286928210 US 



tube and Q is the volumetric flow rate. In fact, given the velocity profile within the tube, 
little flow is actually moving at this speed. 

Turbulent pipe flows 12 are highly complex flows. Predicting the details of any 
turbulent flow is problematic, however, much is known regarding the statistical properties 
5 of the flow. For instance, turbulent flows contain self-generating, coherent vortical 

structures often termed "turbulent eddies". The maximum length scale of these eddies is set 
by the diameter of the pipe 14. These structures remain coherent for several tube diameters 
downstream, eventually breaking down into progressively smaller eddies until the energy is 
dissipated by viscous effects. 

10 Experimental investigations have established that eddies generated within turbulent 

boundary layers convect at roughly 80% of maximum flow velocity. For tube flows, this 
implies that turbulent eddies will convect at approximately the volumetrically averaged 
flow velocity within the pipe 14. The precise relationship between the convection speed of 
turbulent eddies and the flow rate for each class of meters can be calibrated empirically as 

1 5 described below. 

The apparatus 170 of Fig. 15 determines the convection velocity of the vortical 
disturbances within the flow by cross correlating imsteady pressure variations using an 
array of unsteady pressure sensors, similar to that shown in U.S. Patent Application Serial 
No. 10/007,736, filed November 8, 2001, entitled "Flow Rate Measurement Using Unsteady 

20 Pressures", which is incorporated herein by reference. 

Referring to Fig. 15, the apparatus 170 includes a sensing section 172 along a pipe 
14 and a signal processing unit 174. The pipe 14 has two measurement regions 176,178 
located a distance AX apart along the pipe 14. At the first measurement region 176 are two 
unsteady (or dynamic or ac) pressure sensors 180,182, located a distance Xi apart, capable 

25 of measuring the xmsteady pressure in the pipe 14, and at the second measurement region 
178, are two other unsteady pressure sensors 84,86, located a distance X2 apart, capable of 
measuring the unsteady pressure in the pipe 14. Each pair of pressure sensors 180,182 and 
184,186 act as spatial filters to remove certain acoustic signals from the unsteady pressure 
signals, and the distances X|,X2 are determined by the desired filtering characteristic for 

30 each spatial filter, as discussed hereinafter. 



CiDRA Corporation Proprietary 
CC-0676 



-22- 



EV 286928210 US 



The apparatus 170 of the present invention measures velocities associated with 
unsteady flow fields and/or pressure disturbances represented by 188 associated therewith 
relating to turbulent eddies (or vortical flow fields), inhomogeneities in the flow, or any 
other properties of the flow, liquid, vapor, or pressure, having time varying or stochastic 
5 properties that are manifested at least in part in the form of unsteady pressures. The vortical 
flow fields are generated within the flow of the pipe 14 by a variety of non-discrete sources 
such as remote machinery, pumps, valves, elbows, as well as the fluid or mixture flow itself. 
It is this last source, the fluid flowing within the pipe, that is a generic source of vortical 
flow fields primarily caused by the shear forces between the flow 12 and the wall of the 

10 tube that assures a minimum level of disturbances for which the present invention takes 
unique advantage. The flow generated vortical flow fields generally increase with mean 
flow velocity and do not occur at any predeterminable fi-equency. As such, no external 
discrete vortex-generating source is required within the present invention and thus may 
operate using passive detection. It is within the scope of the present that the pressure sensor 

15 spacing may be known or arbitrary and that as few as two sensors are required if certain 
information is known about the acoustic properties of the system as will be more fully 
described herein below. 

The vortical flow fields 188 are, in general, comprised of pressure disturbances 
having a wide variation in length scales and which have a variety of coherence length scales 

20 such as that described in the reference "Sound and Sources of Sound", A. P.Dowling et al, 
Halsted Press, 1983, which is incorporated by reference to the extend of understanding the 
invention. Certain of these vortical flow fields 188 convect at or near, or related to the 
mean velocity of at least one of the elements within a mixture flowing through the pipe 14. 
The vortical pressure disturbances 188 that contain information regarding convection 

25 velocity have temporal and spatial length scales as well as coherence length scales that 

differ from other disturbances in the flow. The present invention utilizes these properties to 
preferentially select disturbances of a desired axial length scale and coherence length scale 
as will be more fully described hereinafter. For illustrative purposes, the terms vortical flow 
field and vortical pressure field will be used to describe the above-described group of 

30 unsteady pressure fields having temporal and spatial length and coherence scales described 
herein. 
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Also, some or all of the functions within the signal processing unit 174 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 
implemented using analog and/or digital hardware, having sufficient memory, interfaces, 
and capacity to perform the functions described herein. 
5 In particular, in the processing unit 174, the pressure signal Pi(t) on the line 190 is 

provided to a positive input of a summer 200 and the pressure signal P2(t) on the line 191 is 
provided to a negative input of the summer 200. The output of the summer 200 is provided 
to line 204 indicative of the difference between the two pressure signals Pi,P2 (e.g., Pi- 

P2=Pasl). 

10 The pressure sensors 180,182 together with the sxmmier 200 create a spatial filter 

176. The line 204 is fed to bandpass filter 208, which passes a predetermined passband of 
frequencies and attenuates frequencies outside the passband. In accordance with the present 
invention, the passband of the filter 208 is set to filter out (or attenuate) the dc portion and 
the high frequency portion of the input signals and to pass the frequencies therebetween. 

15 Other passbands may be used in other embodiments, if desired. Passband filter 208 
provides a filtered signal Pasfl on a line 212 to Cross-Correlation Logic 216, described 
hereinafter. 

The pressure signal P3(t) on the line 192 is provided to a positive input of a summer 
202 and the pressure signal P4(t) on the line 193 is provided to a negative input of the 

20 summer 202. The pressure sensors 83,84 together with the summer 202 create a spatial 
filter 1 78. The output of the simmier 202 is provided on a line 206 indicative of the 
difference between the two pressure signals P3,P4 (e.g., P3-P4= Pas2). The line 206 is fed to a 
bandpass filter 210, similar to the bandpass filter 108 discussed hereinbefore, which passes 
frequencies within the passband and attenuates frequencies outside the passband. The filter 

25 210 provides a filtered signal Pasf2 on a line 214 to the Cross-Correlation Logic 216. The 
signs on the summers 200,202 may be swapped if desired, provided the signs of both 
summers are swapped together. In addition, the pressure signals Pi,P2,P3,P4 may be scaled 
prior to presentation to the summers 200,202. 

The Cross-Correlation Logic 216 calculates a known time domain cross-correlation 

30 between the signals Pasn and Pasc on the lines 212,214, respectively, and provides an output 
signal on a line 218 indicative of the time delay x it takes for an vortical flow field 188 (or 
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vortex, stochastic, or vortical structure, field, disturbance or perturbation within the flow) to 
propagate from one sensing region 176 to the other sensing region 178. Such vortical flow 
disturbances, as is known, are coherent dynamic conditions that can occur in the flow which 
substantially decay (by a predetermined amount) over a predetermined distance (or 
5 coherence length) and convect (or flow) at or near the average velocity of the fluid flow. As 
described above, the vortical flow field 188 also has a stochastic or vortical pressure 
disturbance associated with it. In general, the vortical flow disturbances 188 are distributed 
throughout the flow, particularly in high shear regions, such as boundary layers (e.g., along 
the inner wall of the tube 14) and are shown herein as discrete vortical flow fields 188. 

10 Because the vortical flow fields (and the associated pressure disturbance) convect at or near 
the mean flow velocity, the propagation time delay t is related to the velocity of the flow by 
the distance AX between the measurement regions 176,178, as discussed hereinafter. 

Referring to Fig. 15, a spacing signal AX on a line 220 indicative of the distance AX 
between the sensmg regions 176,178 is divided by the time delay signal x on the line 218 by 

15 a divider 222 which provides an output signal on the line 196 indicative of the convection 
velocity Uc(t) of the saturated vapor/liquid mixture flowing in the pipe 14, which is related 
to (or proportional to or approximately equal to) the average (or mean) flow velocity Uf(t) 
of the flow 12, as defined below: 

Uc(t) = AX/t oc Uf(t) Eq. 1 

20 

The present invention uses temporal and spatial filtering to precondition the pressure 
signals to effectively filter out the acoustic pressure disturbances Pacoustic and other long 
wavelength (compared to the sensor spacing) pressure disturbances in the tube 14 at the two 
sensing regions 176,178 and retain a substantial portion of the vortical pressure disturbances 

25 Pvorticai associated with the vortical flow field 188 and any other short wavelength (compared 
to the sensor spacing) low frequency pressure disturbances Pother- In accordance with the 
present invention, if the low fi-equency pressure disturbances Pother are small, they will not 
substantially impair the measurement accuracy of Pvorticai- 

The second technique of determining the convection velocity of the vortical 

30 disturbances within the flow 12 is by characterizing the convective ridge of the vortical 
disturbances using an array of imsteady pressure sensors, similar to that shown in U.S. 
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Patent Application Serial No. 09/729,994, filed December 4, 2000, entitled "Method and 
Apparatus for Determining the Flow Velocity Within a Pipe", which is incorporated herein 
by reference. 

The sonar flow metering methodology uses the convection velocity of coherent 
5 structure with turbulent pipe flows 12 to determine the volumetric flow rate. The 

convection velocity of these eddies 188 is determined by applying sonar arraying processing 
techniques to determine the speed at which the eddies convect past an axial array of 
unsteady pressure measurements distributed along the pipe 14. 

The sonar-based algorithms determine the speed of the eddies 188 by characterizing 
10 both the temporal and spatially frequency characteristics of the flow field. For a train of 
coherent eddies convecting past a fixed array of sensors, the temporal and spatial fi-equency 
content of pressure fluctuations are related through the following relationship: 



G> 

^convect 

15 Here k is the wave number, defined as k=27c/A- and has imits of 1 /length, © is the temporal 
fi-equency in rad/sec, and Uconvect is the convection velocity. Thus, the shorter the 
wavelength (larger k) is, the higher the temporal firequency. 

In sonar array processing, the spatial / temporal fi-equency content of time stationary 
soimd fields are often displayed using "k-co plots". K-co plots are essentially three- 

20 dimensional power spectra in which the power of a sound field is decomposed into bins 
corresponding to specific spatial wave numbers and temporal firequencies. On a k-© plot, 
the power associated with a pressure field convecting with the flow is distributed in regions, 
which satisfies the dispersion relationship developed above. This region is termed **the 
convective ridge" (Beranek, 1992) and the slope of this ridge on a k-w plot indicates the 

25 convective velocity of the pressure field. This suggests that the convective velocity of 

turbulent eddies, and hence flow rate within a tube, can be determined by constructing a k-co 
plot from the output of a phased array of sensor and identifying the slope of the convective 
ridge. 

Fig. 16 shows an example of a k-co plot generated firom a phased array of pressure 
30 sensors. The power contours show a well-defined convective ridge. A parametric 
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optimization method was used to determine the "best" Hne representing the slope of the 
convective ridge 200. For this case, a slope of 14.2 ft/sec was determined. The 
intermediate result of the optimization procedure is displayed in the insert, showing that 
optimized value is a unique and well-defined optima. 
5 The k-w plot shown in Fig. 16 illustrates the fundamental principle behind sonar 

based flow measure, namely that axial arrays of pressure sensors can be used in conjunction 
with sonar processing techniques to determine the speed at which naturally occurring 
turbulent eddies convect within a pipe. 

The present invention will now be described with reference to Fig. 17 wherein the 

10 discussions based on the calculation of various parameters and properties are detailed herein 
above with reference to the various Figures. In accordance with the present invention 
utilizing a probe 10,1 10,170 to determine the speed of soxmd propagating through a flow 12, 
such as a mixture, provides various specific properties of a mixture and the velocity of the 
mixture and further utilizing logic comprising information about the flow 12 based on the 

15 measured parameters. The steady state pressure and temperature of the mixture may be 

measured by any known or contemplated method as represented by 270 from which various 
fluid properties may be determined from tables or graphs of the known relationships for 
speed of sound and density for the two phases of the mixture as represented by 271 . The 
speed of soimd propagating through the mixture is determined by the apparatus 10,1 10,170 

20 of the present invention as set forth herein above and represented by 272. The quality of the 
saturated vapor/liquid mixture is determined from the fluid properties of 271 combined with 
the mixture speed of soimd 272 using the Wood equation (or similar) as set forth herein 
above and represented by 273. The present invention also enables the determination of 
other properties of the mixture such as enthalpy and density as set forth by 274 by 

25 combining the fluid properties of 27 1 with the quality or composition of the mixture from 
273. The present invention further enables the determination of the velocity of the mixture 
by the methods described herein above as represented by 275. The total volumetric flow 
rate of the mixture is thereby determined as represented by 276 and when combined with 
the parameters of other properties of the mixture such as enthalpy and density as set forth by 

30 274 various flux rates of the mixture such as enthalpy and mass flow rates are enabled as 
represented by 277. 
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Another embodiment of the present invention include a pressure sensor such as pipe 
strain sensors, accelerometers, velocity sensors or displacement sensors, discussed 
hereinafter, that are mounted onto a strap to enable the pressure sensor to be clamped onto 
the pipe. The sensors may be removable or permanently attached via known mechanical 
5 techniques such as mechanical fastener, spring loaded, clamped, clam shell arrangement, 
strapping or other equivalents. These certain types of pressure sensors, it may be desirable 
for the pipe 12 to exhibit a certain amount of pipe compliance. 

Instead of single point pressure sensors 18-21, at the axial locations along the pipe 
12, two or more pressure sensors may be used around the circumference of the pipe 12 at 

10 each of the axial locations. The signals from the pressure sensors around the circumference 
at a given axial location may be averaged to provide a cross-sectional (or circumference) 
averaged unsteady acoustic pressure measurement. Other numbers of acoustic pressure 
sensors and annular spacing may be used Averaging multiple annular pressure sensors 
reduces noises from disturbances and pipe vibrations and other sources of noise not related 

15 to the one-dimensional acoustic pressure waves in the pipe 12, thereby creating a spatial 
array of pressure sensors to help characterize the one-dimensional sound field within the 
pipe 12. 

The pressure sensors 18-21 of Fig. la described herein may be any type of pressure 
sensor, capable of measuring the unsteady (or ac or dynamic ) pressures within a pipe 14, 

20 such as piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers 
(or geophones), velocity measuring devices, displacement measuring devices, etc. If optical 
pressure sensors are used, the sensors 18-21 may be Bragg grating based pressure sensors, 
such as that described in US Patent Application, Serial No. 08/925,598, entitled " High 
Sensitivity Fiber Optic Pressure Sensor For Use In Harsh Environments", filed Sept. 8, 

25 1997, now U.S. Patent 6,016,702, and in US Patent Application, Serial No. 10/224,821, 
entitled " Non-Intrusive Fiber Optic Pressure Sensor for Measuring Unsteady Pressures 
within a Pipe", which are incorporated herein by reference. In an embodiment of the 
present invention that utilizes fiber optics as the pressure sensors 14 they may be connected 
individually or may be multiplexed along one or more optical fibers using wavelength 

30 division multiplexing (WDM), time division multiplexing (TDM), or any other optical 
multiplexing techniques. 
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It is also within the scope of the present invention that any other strain sensing 
technique may be used to measure the variations in strain in the tube, such as highly 
sensitive piezoelectric, electronic or electric, strain gages attached to or embedded in the 
tube 14. 

5 In certain embodiments of the present invention, a piezo-electronic pressure 

transducer may be used as one or more of the pressure sensors 15-18 and it may measure 
the unsteady (or dynamic or ac) pressure variations inside the tube 14 by measuring the 
pressure levels inside of the tube. In an embodiment of the present invention, the sensors 14 
comprise pressure sensors manufactured by PCB Piezotronics. In one pressure sensor there 

10 are integrated circuit piezoelectric voltage mode-type sensors that feature built-in 

microelectronic amplifiers, and convert the high-impedance charge into a low-impedance 
voltage output. Specifically, a Model 106B manufactured by PCB Piezotronics is used 
which is a high sensitivity, acceleration compensated integrated circuit piezoelectric quartz 
pressure sensor suitable for measuring low pressure acoustic phenomena in hydraulic and 

15 pneumatic systems. It has the unique capability to measure small pressure changes of less 
than 0.001 psi under high static conditions. The 106B has a 300 mV/psi sensitivity and a 
resolution of 91 dB (0.0001 psi). 

The pressure sensors incorporate a built-in MOSFET microelectronic amplifier to 
convert the high-impedance charge output into a low-impedance voltage signal. The sensor 

20 is powered fi-om a constant-current source and can operate over long coaxial or ribbon cable 
without signal degradation. The low-impedance voltage signal is not affected by 
triboelectric cable noise or insulation resistance-degrading contaminants. Power to operate 
integrated circuit piezoelectric sensors generally takes the form of a low-cost, 24 to 27 
VDC, 2 to 20 mA constant-current supply. A data acquisition system of the present 

25 invention may incorporate constant-current power for directly powering integrated circuit 
piezoelectric sensors. 

Most piezoelectric pressure sensors are constructed with either compression mode 
quartz crystals preloaded in a rigid housing, or unconstrained tourmaline crystals. These 
designs give the sensors microsecond response times and resonant fi-equencies in the 

30 hundreds of kHz, with minimal overshoot or ringing. Small diaphragm diameters ensure 
spatial resolution of narrow shock waves. 
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The output characteristic of piezoelectric pressure sensor systems is that of an AC- 
coupled system, where repetitive signals decay until there is an equal area above and below 
the original base line. As magnitude levels of the monitored event fluctuate, the output 
remains stabilized aroimd the base line with the positive and negative areas of the curve 
5 remaining equal. 

It is also within the scope of the present invention that any strain sensing technique 
may be used to measure the variations in strain in the pipe, such as highly sensitive 
piezoelectric, electronic or electric, strain gages and piezo-resistive strain gages attached to 
the pipe 12. Other strain gages include resistive foil type gages having a race track 

10 configuration similar to that disclosed U.S. Patent Application Serial No. 09/344,094, filed 
June 25, 1999, now US 6,354,147, which is incorporated herein by reference. The 
invention also contemplates strain gages being disposed about a predetermined portion of 
the circumference of pipe 12. The axial placement of and separation distance AXi, AX2 
between the strain sensors are determmed as described herein above. 

15 It should be imderstood that any of the features, characteristics, alternatives or 

modifications described regarding a particular embodiment herein may also be applied, 
used, or incorporated with any other embodiment described herein. 

Although the invention has been described and illustrated with respect to exemplary 
embodiments thereof, the foregoing and various other additions and omissions may be made 

20 therein and thereto without departing from the spirit and scope of the present invention. 
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